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We analyze pseudogap phenomena widely observed in the underdoped cuprates. We assume 
the existence of a strong d-wave pairing force competing with antiferromagnetic(AFM) fluctua- 
tions and the formation of flat and damped dispersion around the (it, 0) and (0, 7r) region as two 
important elements caused by the proximity from the Mott insulator. Using the mode-mode 
coupling theory for the d-wave superconducting(dSC) and AFM fluctuations, we reproduce 
basic properties of the pseudogap seen in the magnetic resonance, neutron scattering, angle 
resolved photoemission and tunneling measurements in the cuprates. Then minimal require- 
ments to understand the pseudogap phenomena are clarified as the above two elements. A 
strong competition of the pairing with the antiferromagnetic fluctuations suppresses the transi- 
tion temperature thereby generates the pseudogap in the underdoped region while the weakness 
of the AFM fluctuations leads to the absence of the pseudogap at the optimal doping concen- 
tration. 

KEYWORDS: high-T c superconductivity, pseudogap, antiferromagnetism, d x 2_ y 2-wa\/e superconductivity, Mott 
transition, strong-coupling superconductor, mode-mode coupling 



§1. Introduction 

Pseudogap is one of the most remarkable phenomenon in the underdoped region of high-T c 
cuprates. It is observed both in spin and charge excitations in which gap structure emerges from a 
temperature Tpg well above the superconducting transition point T c . The gap structure is observed 
in various different experimental probes such as NMR relaxation time, the Knight shift, neutron 
scattering, tunnenling, photoemission, specific heat, optical conductivity, and DC resistivity .B* 
The angle resolved photoemission spectra (ARPES)ii have revealed that the pseudogap starts 
growing first in the region around (tt,0) and (0, it) from T = Tpq much higher than T c . In the 
earlier work^ these momentum regions are known as the region where the quasiparticle dispersion 
becomes unusually fiat and strongly damped. This gap structure continuously merges into the 
d x 2_ y 2 gap below T c . We call such (it, 0) and (0, it) momenta "flat spots" and the region around them 
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"flat shoal region" . This region is also known to be particularly important in the understanding of 
the metal-insulator transition and its scaling properties!^ 

One puzzling experimental observation is that the pseudogap structure appears in 
l/TiT^'i'i'00 while in many cases I/T2G continuously increases with the decrease in temper- 
ature with no indication of the pseudogap. In addition, the so called resonance peak appears in 
the neutron scattering experiments. A resonance peak sharply grows at a finite frequency below 
T c with some indications even at T c < T < Tpq. This peak frequency iv* decreases with lowering 
doping concentration implying a direct and continuous evolution into the AFM Bragg peak in the 
undoped compounds. The neutron and T2G data support the idea that the AFM fluctuations are 
suppressed around u = but transferred to a nonzero frequency below Tpq- These observations 
require the framework treating the superconducting and AFM fluctuations on an equal footing. 

As we discuss in this paper, the pseudogap phenomena are well understood as a consequence of 
two fundamental aspects of the cuprates. The first is proximity effects from the Mott insulator 
near the metal-insulator transition where strong Coulomb repulsion generates a strong and critical 
momentum and energy dependences in electron excitations. The coherence temperature (effective 
Fermi energy) is unusually suppressed due to this proximity. In the momentum space, the flat shoal 
region appears and this region determines the basic character of the metal-insulator transition. 
The strong correlation effects appear most critically in this region with formation of flattened and 
strongly damped dispersion. Because of its flatness with diverging density of states, the doping 
effects are determined predominantly by this region. The flat shoal region has a fundamental 
importance also in clarifying the mechanism of the pseudogap since the ARPES result shows that 
the pseudogap is first formed from this region. 

The second fundamental aspect in the cuprates is the strong coupling nature of the pairing in- 
teraction. The short coherence length observed in the cuprate superconductors implies that the 
effective Fermi energy Ep is comparable or even smaller than the energy scale of the pairing interac- 
tion in contrast with the conventional BCS superconductors. This is indeed a natural consequence 
of the suppressed coherence and Ep for a metal near the Mott insulator when the pairing force is 
kept constant. The appearance of the pseudogap region characterized by a separation of supercon- 
ducting transition temperature T c and the onset of pairing fluctuation is naturally understood from 
this strong coupling character. In this paper, we see that this separation is strongly enhanced by the 
repulsive mode-mode coupling between dSC and AFM fluctuations because the AFM fluctuations 
suppress superconducting T c . 

To reach full understanding of the pseudogap phenomena, we need satisfactory descriptions of 
both of the above two aspects, although even a complete description of each aspect alone has never 
been given in the literature. The formation of the flat shoal region is observed in numerical stud- 
iesl'lliP while it is not well reproduced in self-consistent treatments by the diagrammatic approaches. 
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Concerning the second aspect, several strong coupling approaches for the superconducting phase 
in the literature have not seriously treated neither the competition of the pairing interaction with 
the antiferromagnetic fluctuations nor the suppression of the coherence temperature. These so far 
ignored elements are actually crucial elements in realizing a region of strong pairing fluctuation 
above the suppressed T c . Furthermore, under the serious competition of the two fluctuations, the 
origin of the strong pairing force has not been clarified yet. This is presumably because the pair- 
ing force itself results from the incoherent high-energy excitations under the proximity effect from 
the Mott insulator while microscopic theory for such incoherent part is not fully developed. This 
difficulty indeed becomes clear when we see the results of this paper, where the antiferromagnetic 
fluctuations must be repulsive with the pairing force at low energies in reproducing the pseudogap 
formation and a rather high-energy incoherent excitations are required for the origin of the pairing 
force. 

Keeping the present stage of the above understanding in mind, we develop a theoretical framework 
of the pseudogap phenomena to account for all of the basic experimental results. The scope of this 
paper is not so ambitious that a full microscopic theory is constructed. In this paper we rather aim 
at giving a framework where the various experimental results are reproduced from a theory starting 
from the observation of the above two aspects, namely the strong coupling nature of pairing and 
critically strong momentum dependence in the quasi-particle excitations. In this paper, irreducible 
pairing interaction is rather given as input and assumed to exist. We also pick up enhanced 
contributions from the flat shoal region to the dSC and AFM susceptibilities and impose a cutoff 
in the integral over momentum space to mimic the dominant contribution of this region, although 
the flatness and the damping are not completely expressed in this scheme. We then construct a 
mode-mode coupling theory for both the d-wave pairing and antiferromagnetic fluctuations on an 
equal footing. We will show that such minimal requirements are enough to reproduce the basic 
experimental results of the pseudogap phenomena.0El) 

Although the existence of the flat shoal region plays main role for the criticality of the metal- 
insulator transition and the formation of the pseudogap, a subtlety arises in some physical quantities 
for the role of the other region as (7r/2,7r/2) point. In fact, the DC transport and damping of 
the magnetic excitations could substantially be influenced from doped holes in the other dispersive 
region. Although the transition to the Mott insulator is not accompanied by the critical behavior of 
the relaxation time r but by r-independent quantities as the Drude weight and the compressibility, 
the noncritical quantities such as the DC transport and the magnetic relaxation may sensitively 
depend on r. This is particularly true for the damping of the magnetic excitations under the 
pseudogap formation. If contributions from the (-7r/2,7r/2) region would be absent, the damping of 
the magnetic excitation would be strongly reduced when the pseudogap is formed around (vr,0). 
However, under the pseudogap formation, the damping can be determined by the Stoner continuum 
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generated from the (jr/2, n/2) region and can remain constant. This process is in fact important if 
the quasiparticle damping around the (tt/2,tt/2) region is large as in the case of La 214 compounds. 
Since the whole momentum dependence of the quasiparticle damping is not easy to derive in the 
present stage, and the damping of the magnetic excitations are determined from rather complicated 
combination from the both flat and dispersive regions, in this paper, we leave the damping of the 
magnetic excitations as an input from the outside of the framework based on phenomenological 
grounds. The formation of the pseudogap itself is a rather universal consequence of the strong 
coupling superconductors. However, as we see below, the actual behavior may depend on this 
damping. For example, we show below that the damping generated by the (tt/2,tt/2) region 
sensitively destroy the resonance peak structure observed in the neutron experimental results. 

§2. Mode-mode coupling treatment 

We consider a 2D strongly correlated electron system and treat AFM and pairing fluctuations 
simultaneously. 0) Following the argument in §1, we represent the partition function of the system 
by the functional integral over both of the AFM and dSC auxiliary fields, 4> a and </>d introduced 
by the Stratonovich-Hubbard transformation. 

After integrating out the fermions degrees of freedom, the following effective action is obtained, 

S = S® + + S® + sW + fig) + s% , (2.1) 

s l 2) =/ ? E/ dUx^K.gj^K.Ql'if-i^,-?), (2-2) 

n J 

S d ] = 0^2 d 2 q^ d xl 1 (iwn,q)0d(i^n,q)0d(iwn,q), (2.3) 

n •* 

S$ = 0u aa J2 / d2< ?i / d 2 9 2 / ^Qa^A^n^qi)-^^^, q 2 )4>A^n :i , q 3 )-4>a(^n 4 , q^), (2.4) 

ni,n 2 ,n 3 J J J 

S dd = P u ddYl / d2 9l / d 2 <? 2 / d 2 g3^d(-iWni,-qi)0d(iWn 2 ,q2)^d(-iWn3,-q 3 )0d(iWn4,94X2.5) 

S ad = 2 P u adJ2 / d 2 9i / d 2 Q 2 / d 2 Q 3 CT (i^n 1 ,qi)-^ (J (iw„ 2 ,q2)0d(-i^n3,-q3) ( ? !> d(iw„ 4 ,q4),(2.6) 

where 4> a is the three-component vector field corresponding to the spin, and 4>d and (f>d are the 
pairing fields creating and annihilating a pair of electrons, respectively. (3 is the inverse temperature 
and u n and Q m are bosonic and fermionic Matsubara frequencies. Here 

X^ n , q )=Aj^ + ( q -Qf+^ + 4), 1 (2-7) 

C \ A ( t(0)-2 i 2 , 7dkn| . Wn\ , s 

Xd(iuJn,q)=A d [C d J +q + — 2— + — ' ( 2 - 8 ) 
V c d c d J 
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are AFM and dSC dynamical susceptibilities, q 4 = — q 1 — q 2 — q%, = —n\ — n 2 — n 3 , and Q 
is the AFM ordering wave vector. The bare correlation lengths ffi and £^ are determined from 
the inverse of the bare susceptibilities calculated from the bare dispersion e^. Here, we take the 
Hubbard model with nearest-neighbor transfer t and the next-nearest-neighbor transfer t' , which 
leads to = —2t(cosk x + cosk y ) — At' (cosk x cosk y + 1) — //. The chemical potential \i is measured 
from the flat spots. We take the form 

e ] - 2 « i - ^ log — ^-ty iog f n ^ 

t max{fi,t,l\ max{fi,l\ 



Vt 2 - W 2 T max{/j,, T\ 



er « i - iog ^ bg - (2.io) 



which is valid for the contributions from the (7r, 0) and (0, tt) regions. Here E c is a ultraviolet cutoff 
in the energy scale of the bandwidth and we have included coefficients of the double logarithms 
into the original AFM and dSC Gaussian coupling coefficients F a and i~d- In the above forms, the 
role of the flat shoal region is emphasized and selectively picked up by a cutoff imposed in the 
momentum space. 

We note that although we calculate the bare correlation length from the Hubbard model, our 
scheme is beyond the scope of it because we introduced the AFM and dSC coupling constants u, 
r, A a and Ad as phenomenological parameters. 

In this paper we restrict ourselves to the AFM ordering vector at the commensurate value (tt, tt). 
We also neglect possible long-range features of Coulomb interaction which may lead to gapful 
dSC excitations instead of the Goldstone mode even in the dSC ordered state, as in the s-wave SC 
state. We have confirmed that the susceptibility Xd taken to satisfy the Anderson-Higgs mechanism 
instead of (2.8) dose not alter the qualitative feature of the pseudogap behavior obtained in this 
paper. The phase excitations (Higgs bosons) are not treated separately from the amplitude modes. 
The velocity of spin and pairing collective modes are denoted by c a and Cd • The damping constants 
are given by 7d and j a . 

Following the argument in §1 for the damping of the AFM fluctuations, we introduce a phe- 
nomenological form of 7o- and 7d- The origin of 7 is mainly from continuum of the Stoner excita- 
tions and the amplitude strongly depends on low-energy quasiparticle excitations. This low-energy 
part of damping becomes negligible if some kind of long-ranged order appears. It may also be sup- 
pressed if the correlation length gets longer. When only one type of fluctuations with the correlation 
length £ exists, a plausible dependence for long £ would be 7 = 7 < ° ) £~ v '. In case of the enhanced 
d-wave correlation length, however, the situation is not so simple, because the low-energy excita- 
tions around (vr/2, tt/2) are not suppressed due to the nodal structure of the d-wave gap. To get 
qualitative results, these suggest us a rough form for the damping as 7 oc 7i/(£ ^ + +72/£ct ¥ '- 
Here the first term represents the contribution from the (tt, 0) region and the second term is from 
the (tt/2, tt/2) region. Since the damping does not have critical change in the pseudogap region, 
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the above rough form may be a good starting point to get an idea for its role. Depending on the 
relative amplitude of the first and second terms, we may take only the dominant term between 
these two. In the pairing dominant region as in the subject of this paper, these two choices may 
also be expressed simply as 

7^ = 27^/^ + ^), (2-11) 

where we should take tp = in the latter choice for the (tt/2, tt/2) contribution. In terms of bosonic 
excitations, the relaxation times of collective modes should be determined by the time necessary 
to propagate the scale of the longest correlation length, because the damping is not effective as 
far as the excitations are propagating inside such an ordered domain. Thus we take tp = 1 for the 
former choice. The case tp = 1 represents the one where the damping at (tt, 0) is overwhelming over 
the (tt/2, 7r/2) region thus is generally adequate for the underdoped region. The exception is the 
La-based 214 compounds, where the quasiparticle damping around (tt/2, tt/2) is unusually large 
presumably because of charge ordering fluctuations. The optimally doped compounds are rather 
expressed by tp = because the relatively weak damping around (tt, 0) does not allow us to neglect 
the contribution from the (tt/2, tt/2) region any more. 

Although we take different choices, (p = or 1 for optimal and underdoped regions, respectively, 
we note that the form for the damping does not alter the formation of the pseudogap itself. The 
pseudogap formation is a consequence of large d-wave coupling constant competing with AFM 
fluctuations under unusually suppressed coherence temperature. We will see later that the pseu- 
dogap appears only in the underdoped situation in our scheme. This is simply because the mutual 
competition between dSC and AFM fluctuations are severe there and additionally the coherence 
temperature becomes low. The damping form is crucial only for the appearance of the resonant 
peak in the underdoped region. If we take ip = instead of our present choice tp = 1 in the 
underdoped region, the pseudogap survives but the resonant peak would not appear. 

Using the effective action obtained above, we perform renormalization process for the mode-mode 
coupling terms. This is a similar procedure to the SCR theory developed by Moriya and coworkers 
for spin fluctuations© In our case the mode-mode coupling terms consist of those between AFM 
and AFM fluctuations, dSC and dSC fluctuations, and AFM and dSC fluctuations. Following the 
mode-mode coupling scheme, £ CT and £d are determined selfconsistently from 



= d 0) ~ 2 + J o —J o coth [*wImx CT (w, k) + u CTd Imxd(w, k)] , (2.12) 

£d 2 = Q + / — / TT^f coth krdlmx^w, k) + Wmxd(w, *)] , (2-13) 
Jo tt Jo \2iry 11 

where in the susceptibilities \: the bare correlation lengths £o- and in (2.7) and (2.8) are 
replaced with the renormalized ones without (0). 
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In this formalism for two-dimensional systems, the system can be ordered only at T = 0. It 
agrees with the Mermin- Wagner theorem® for the AFM order. For the dSC, however, the K-T 
transition at nonzero temperatures is not reproduced since the SCR theory cannot describe such 
topological K-T transition at Tkt- However, the renormalized superfluid stiffness determines the 
temperature scale where the pairing correlation length starts growing strongly. Since Tkt is of the 
order of the stiffness according to the K-T theory Tkt is close to this crossover temperature T* 
in our theory, below which the spin correlation length starts decreasing. In our analysis, we take 
this temperature scale as the signature of the K-T transition. 

For our calculation, we choose two sets of parameter values. In one, parameter values for typical 
underdoped compounds such as YBa2Cus06.63 and YBa2Cu40s are taken and for the other, typical 
optimally doped case such as YBa2Cus07 is considered. For details of the determination of the 
parameters readers are referred to the paperBH) 

§3. Results 

First, we discuss results for optimally doped systems. The parameter values are chosen from 
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Fig. 1. The spin correlation length (solid line), 1/TiT (dashed line) normalized by its value at T = t, and dSC 
correlation length (long-dashed line) plotted as functions of temperature. The parameter values are for the optimally 
doped cuprates and given as \F a \ — |7d| = 0.15, u ac , — 2.6, u CT d = 1-2, u^d = 1-8, /i = 0.03 and t' = —0.03. This 
choice corresponds to YBa2Cu 3 07- 



an optimally doped compound, YBa2Cu307LP and (f = is also taken. We plot the calculated 
dSC correlation length £<j, 1/T\T and spin correlation length £ CT in Fig. 1. In this case, the spin 
correlation length first increases down to the temperature T* and then decreases with further 
decrease in temperature. It makes a crossover between the regime T < T* dominated by the dSC 
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renormalized classical fluctuations and the thermally fluctuating regime T > T*. 

The spin-lattice relaxation rate 1/Ti of a 63 Cu nuclei is evaluated from the momentum sum of the 
imaginary part of the dynamical spin susceptibility. Here we have not considered the k-dependence 
of the nuclear form factor seriously, because it does not alter the basic feature. We see that \jT\T 
has basically the same temperature dependence as Our present result for (p = is also similar 
to the results for Tl2Ba2Cu06, or HgE^CuC^+tfS Our results are totally consistent with the 
absence of the pseudogap region seen in experimental results of optimally as well as overdoped 
cuprates. 

We next consider the underdoped region with a special emphasis on the resonance peak behavior. 
The calculated £ CT , ^ and \jT\T are shown in Fig. 2 for the parameter values corresponding to 




Fig. 2. Temperature dependence of the spin correlation length (solid line), \jT\T (dashed line), and the dSC 
correlation length (long-dashed line) in the ip — 1 case. We have taken \r a \ — \r c \\ — 0.7, u aa = 3.14, iidd = 3.05, 
u a d = 1.0, /i = 0.02 and t' = —0.02, which correspond to an underdoped cuprate, YBa2Cu3C>6.63- As in Fig. 1, 
the data of 1/TiT are normalized by its value at T — t. The fact that £d starts increasing rapidly with decrease in 
temperature at T, can be confirmed by the vanishing amplitude of 1/TiT oc £&/(£o + fa) reached near T*. 



an underdoped cuprates, YBa2Cu30g.63- Here we take cp = 1. In contrast to the optimally doped 
case, the spin correlation length has its maximum at a temperature well above T*. With decrease 
in temperature below Tpq, ^ starts growing quicker than This competition between ^ and £ CT 
is an origin of the pseudogap formation. Below T*, the dSC fluctuations go into the renormalized- 
classical regime, which signals the decrease in £ a . We again interpret T* as the rough estimate of 
T c . These properties are also similar to experimental data in underdoped cuprates with a pseudo 
spin gap, such as YBa2Cu40g& and Bi2Sr2CaCu20g.B 
The growth of the pairing correlation also drives reduction of the damping j a in spin excitations 
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and makes underdamped resonance peak at a finite frequency uj = uj* in S(Q,u>). Namely, the 
spectral weight start transferring from u = to the peak region around uj* . A similar crossover 
was previously obtained in a numerical calculation near the quantum transition point between 
dSC and AFM ordered phasesSBil' The peak structure in S(Q,u>) around u>* reproduces some 
qualitative feature in the resonance peak observed experimentallylll'B as we see in Fig. 3. In our 
treatment, uj* is self-consistently determined from the competition between dSC and AFM and the 
value uj* is characterized by the dSC gap amplitude. The AFM fluctuations are pushed out from 
the region lower than uj* due to the triplet excitation gap generated by the dSC gap formation. 




Fig. 3. The spin structure factor S(ui, Q) with Q — (n, n) for the underdoped case with <p = 1. From the overdamped 
side, T/t = 0.069, 0.06, 0.051, 0.042, 0.033 and 0.024. The condition ip = 1 makes the spin excitation underdamped 
and brings about the shift of the low-energy spectral weights to higher energies. It evolves as a resonance peak as 
observed in neutron scattering experimants. 



In our framework of (2.7), uj* has to be proportional to for small j a , while experimentally, 
growth of the correlation length takes place with a fixed finite uj* when the temperature is lowered. 
To reproduce the temperature dependence in experiments, we need to modify the assumed form 
(2.7) as discussed before@0) and will need to consider dominat incoherent part in addition to the 
coherent response. 

We next discuss the single-particle spectral weight© Here we compare our results with ARPES 
data in Bi2212 with similar values for Tpq(~ 170K) and T c (~ 83K) to those in YBa2Cu30g.63- 
We calculate the electronic spectra ImG(w, k) using the same parameter values for the above 
choice of the underdoped cuprates. The single-particle Green's function is defined by G(uj,k) = 
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1/(lo — e(k) — k)). Here we calculate the self-energy within the 1-loop level using 

d 2 fc' rdu' 



l^(u,k)=J^I^-lrnG^,k') 



r*lmxa(u -uj',k- k') (^coth + tanh —J 

r2 g{kf+g(k')\ , 
+ r d Imxd(^ + uj ,k + k) 

( oj + {J uj' 

x coth — — tanh — — 

V 2T 2T 



(3.1) 



with the bare Green's function G^°\u>, k) and g(k) = (cosk x — cosk y )/2. Here C^ -2 has been 
replaced with £~ d - For the prefactor A d in (2.8), we take A d = 4i _1 to give a proper value for the 
midpoint shift in ARPES intensity in the pseudogap region.& 

Figure 4(a) shows ImG(a;, kp) with kp near the flat spot at various temperatures. At the highest 
temperature T > 0.069i, we have a peak at u) = 0, though it is damped by thermal fluctuations. At 
lower temperatures still above Tpq, only the low-energy spectral weights gradually start decreasing. 
We note that dSC correlations grow more rapidly than those of AFM below Tj>q (~ 0. They 
suppress only the low-energy part of the peak in the spectral weights. Well below Tpq, ImG shows 
further loss of weights around u = 0. For the same momentum fcp, we also plot the intensities 
I(uj, kp) = ImG(w, kp)f(uj) in Fig. 4(b), where / is the Fermi function. The energy of the midpoint 
is nearly zero at T = 0.06t(~ Tpq). For T < Tpq, the midpoint shifts to higher binding energies. 
This shift amounts to 0.045i ~ llmeV at T = 0.042t(~ 122K), in agreement with experiments. & 

The momentum dependence shows that the low-energy part of the single-particle excitations 
is under a stronger suppression near the flat spot , while those closer to the nodes are better 
understood as quasiparticles. The calculated results clearly show the formation of the pseudogap 
first from the (n, 0) region. The overall result qualitatively well captures the emergence of the 
pseudogap structure observed in angle-resolved photoemission experiments. 

§4. Summary and Discussion 

By assuming the d-wave attractive channel and the presence of strongly renormalized flat quasi- 
particle dispersion around the (tt, 0) region, we have considered the mode-mode coupling theory 
for the AFM and dSC fluctuations. The pseudogap in the high-T c cuprates is reproduced as the 
region with enhanced dSC correlations and is consistently explained from precursor effects for the 
superconductivity. The existence of the flat shoal region plays a role to suppress the effective Fermi 
temperature Ep. This suppressed Ep and relatively large pairing interaction both drive the 
system to the strong coupling region thereby leads to the pseudogap formation. The pseudogap 
formation is also enhanced by the AFM fluctuations repulsively coupled with dSC fluctuations. 

The pseudogap formation clarified from the interplay of AFM and dSC is summarized as fol- 
lows: When the dSC correlation grows faster but competes severely with the low-energy AFM 
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Fig. 4. Temperature dependence of (a) the imaginary part of the Green's function and (b) the ARPES intensity 
for &f = (7r, 37r/64) for tp = 1. This momentum point is on the Fermi surface and the closest to (n, 0) in our 
calculation. Temperatures in the plotted data are 0.069, 0.06, 0.051, and 0.042 in the energy unit of t from the 
data with larger intensity at co = both for (a) and (b). 

fluctuations, the pseudogap structure appears above T c with a suppression of T c . The pseudogap 
is observed clearly in the suppression of \jT\T. Detailed structure of the pseudogap depends on 
the damping exponent if. By taking a proper choice of parameters for several underdoped cuprates 
with ip = 1, 1/T\T shows a faster decrease at Tpg(> T*) while £ CT continues to increase until T*. 
With this parameter values, the resonance peak at a finite frequency in S(q,u) is also obtained. 
The single-particle spectral weight shows the growth of the gap structure around (it, 0) and (0, n) 
below Tpg- The qualitative similarity between our results for the underdoped case with ip = 1 
and the experimental results in YBa2Cu3C>6.63 and underdoped Bi2Sr2CaCu20§+5 suggests that 
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the damping of the AFM and dSC collective modes decreases in the pseudogap regime at least 
for these compounds. It means that low-energy fermions around the flat spots mainly contribute 
to the damping. This is consistent with the strong damping of quasiparticle around the flat spot 
observed experimentally in the underdoped region. 

If the AFM fluctuations are not strong enough and do not compete severely with the dSC 
fluctuations, the spin correlation length, £<j and l/ 63 TiT both reaches its maximum value only 
at T = T* and then decreases as the temperature decreases, which indicates the absence of the 
pseudogap region. The experimental results in optimally and overdoped cuprates are reproduced 
from this choice of parameters together with </> = 0, namely the case where the damping 7 does 
not depend on the <iSC correlation length. 

The success in reproducing the pseudogap behavior in spin excitations is based on the competition 
between low-energy AFM and dSC fluctuations. Such competition requires the repulsion for u a d > 
0. Then the d-wave attraction cannot be mediated by low-energy spin fluctuations. Although it does 
not necessarily exclude the attraction generated from the high-energy part of the spin fluctuations, 
it requires a formalism for such incoherent contributions beyond the conventional weak-coupling 
approach. 

Further studies are required for a more complete understanding of the pseudogap in the high-T c 
cuprates. Microscopic derivation of our two starting points is the most intriguing future subject. 
We have concentrated on the single-particle excitations only around the flat spots, (tt, 0) and 
(0, 7r). However, in the one-loop level, the origin of the flat dispersion and strong damping in this 
momentum region is not fully clarified. Experimentally the flatness and damping strength appear 
much more pronounced than the expectation from the one-loop analyses. Numerical analyses 
also support that this remarkable momentum dependence around the flat spots is generated by 
the strong correlation effects. We have to calculate self-energy corrections as well as the vertex 
corrections in a self-consistent fashion to clarify the profoundness of such correlation effects. This 
is clearly the step beyond the one-loop level. This will also contribute to clarify how the pairing 
channel appears and how the flat spots are destabilized to the paired singlet. We also note that the 
dominance of the incoherent weight over the quasiparticle weight in the single-particle excitations 
near the metal-insulator transition may require a serious modification in the derivation of the 
AFM and dSC susceptibilities. The Curie- Weiss type form for the dynamic spin susceptibility we 
assumed needs to be reconsidered ,© because the spin susceptibility is also determined mainly from 
the incoherent part of the single-particle excitations which we have not considered at all in this 
paper. 
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